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Abstract 

Recent experimental investigations have revealed a non-Maxwellian 
absorption pattern in the rheological spectra of actin gels, which was 
interpreted in terms of transient bonds. Here we examine the consequences 
of reversible crosslinking on the apparent linear spectra of biopolymer 
solutions theoretically. For a schematic model consisting of a reversibly 
crosslinked power-law fluid we obtain a simple analytical prediction for 
the position of the absorption peak, which is backed up by a numerical 
evaluation of the inelastic glassy wormlike chain model. This establishes 
bond breaking as a nonlinear non-equilibrium effect that can already be 
significant for very small driving amplitudes. Our results may be useful for 
inferring binding affinities and reaction rates of biochemical crosslinkers 
from rheological measurements of in-vitro reconstituted cytoskeletal gels. 

1 Introduction 

Due to rapid recent progress in the biological sciences an increasingly detailed 
understanding of the biochemical response and regulation of cells is quickly 
emerging. On the other hand, also the response of cells to mechanical stimulus 
has been shown to be of crucial physiological importance'^. In this research 
area, sophisticated experimental techniques are at hand'^' and a rich and excit- 
ing phenomenology has been uncovered. Attempts have been made to organize 
the mosaic of observational evidence by identifying universal patterns such as 
power-law rheology'^'ED the pivotal role of prestress in cell mechanics and for 
the physiology of adhering ceUsEM ]^, an( j t ne fluidization response of adhering 
cells upon transient stretchHEMI Despite the advances in the phenomenologi- 
cal characterization, the underlying physical and molecular mechanisms are still 
debatedEEHHl. 

One strategy to unravel the p hysical foundations of cell mechanics is the 
so-called bottom-up approach' 17 * 18 '. It aims at isolating so-called fuctional mod- 
ules of manageable complexity that can be reconstituted and studiedm vitro. 
Examples for the success of this approach range from reconstituted biomimctic 
me mbran e systems^, over cytoskeletal polymer networksl^HII] anc j motility as- 
says' 25 ' 2 ^ to conceptual progress in the understanding of complex in-vivo func- 
tions such as cell adhesion^. 
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It is commonly agreed that the actin cortex is a main contributor to the 
mechanical response of tissue cellsPEH] Therefore, reconstituted actin networks 
with and without crosslinking molecules are studied by many groups and cell- 
like mechanical behavior has indeed been demonstrated'^^Sl An interesting 
feature that has been observed in actin networks crosslinked by a single type of 
crosslinker is the emergence of a distinguished frequency scale in the frequency- 
dependent rheology l 31 * 32 !. A non-Maxwellian relaxation pattern consisting of a 
pronounced peak appears in the loss modulus, accompanied by a shoulder in the 
storage modulus. The increased dissipation has been attributed to the breaking 
of transient bonds and a phenomenological connection between the zero- force off 
rate of the crosslinker and the characteristic frequency seen in rheology has been 
proposed^ 20 * 32 * 33 !. This proposition does not contradict the fact that comparable 
effects are commonly not reported for cells, as in cells there are many different 
types of crosslinkers with a whole distribution of off rates present. 

It therefore seems that a comprehensive theoretical model accounting for 
breaking and reforming of transient bonds could contribute much to the un- 
derstanding of the mechanics of biopolymer networks and of cells. It has been 
argued that a good model should capture the nonequilibrium character of cell 
mechanical processes on top of the polymeric basis of cytoskeletal elasticity 
and bridge the length and time scales between the elementary molecular events 
and meaningful experimental model systems^. This is a precise characteri- 
zation of what the following analytical and numerical analysis aims at. We 
analyze the problem on two levels. First, we study analytically the influence of 
bond-breaking on the nonlinear response of a schematic model for biopolymer 
networks, a transiently crosslinked power-law fluid. We analytically derive an 
absorption pattern in the frequency-dependent rheology and provide a predic- 
tion for its position as a function of parameters such as the equilibrium off rate 
and the relative binding affinity. Second, we report numerical results obtained 
with a recently proposed generic model for the inelastic mechanics of biopoly- 
mer networks, the inelastic glassy wormlike chain^. We find good qualitative 
agreement of both theoretical predictions with each other and with previously 
published results for crosslinked actin networks^ 32 * 33 !. 

2 Models 

2.1 Model for the bond ensemble 

Consider an ensemble of non-interacting bonds which have two well-defined 
states, the bound state at xj, and the unbound state at x u . By convention, 
the energy vanishes in the force-free bound state, and the unbound state has 
energy U. To go from one state to the other, a bond has to pass through a 
metastable transition state of energy £ at x t- External forces / are assumed to 
adiabatically tilt the energy landscap d 36 * 37 l This phenomenological assumption 
can be recovered from more elaborate theories^^ as an approximation limited 
to certain experimental conditions. Beyond that, we do not consider these more 
elaborate theories as particularly helpful for our present task, which is to derive 
a robust practical analytical prediction as to how the rheology of biopolymer 
networks is affected by bond breaking. First, most of them disregard rebinding, 
which is expected to be relevant for us, second they cannot be straightforwardly 
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implemented withour either using a much more elaborate model of the polymer 
network or invoking further speculative assumptions, e.g. on the geometry of 
and the force transduction through the network. We therefore choose to stick 
(for now) to the simple Bell modeP^, in which the rates k + and for binding 
and unbinding are given by 

k + (f) = ke u ~ A ^f, (1) 

and 

k.(f) = ke Ax »f, (2) 

with Ax b / U = \xi,/ u — Xt\, respectively. To simplify the notation, we measure 
energies in units of the thermal energy fcgT, giving lengths and moduli the di- 
mensions of inverse force and and force squared, respectively. The prefactor k 
corresponds to the equilibrium off rate . In the framework of Kramers' theory^, 
it is proportional to the microscopic attempt rate and decreases exponentially 
with the barrier height, k = e _£ T^~ . More sophisticated theoret ical appr oaches 
predict different functional forms of the equilibrium off rate^^t ^ * 41 * 43 * 44 !, also 
including dependencies of the equilibrium off rate on the stiffness of the force 
probe^ 41 * 43 * 44 * and even explici tly take into account a pulling-induced time- 
dependent free-energy landscap d 39 * 40 * 42 * ^. However, at the present state we 
are not interested in the microscopic model details and prefer to treat k as 
an independent parameter, unless we have reason not to do so (further model 
detrails may e.g. matter when one tries to extract biochemical data from rheo- 



logical experiments - see section 3.3 1 



Noting that the fraction of closed bonds and the fraction of open bonds have 
to add to one, the fraction v of bonds in the closed state ("bond fraction") is 
given by the first-order kinetic equation 

v{t) = -{k + + k_)v{t) + k+. (3) 

We are interested in the linear response of the bond fraction to a sinusoidal 
driving force, 

f(t) = fsm(w Q t). (4) 



Expanding the solution of equations ([1|-(|4| with respect to /, truncating at 
linear order, and neglecting transients in the long-time limit, t ^> k(e u + 1), 
yields 

v(t) m v(f = 0) - / ■ {al(uj ) sm(uj t) - a"(u ) cos(w £)), (5) 



with the abbreviations 



, (Ax b + Ax u ) 1 

= " ^ ^2 ( 6 ) 



4cosh^(t//2) - , ( „ \ 

1 + \ k(l+e u ) ) 

and 



„ (Ax b + Ax u ) k(l+eU) 



4cosh^(C//2) , , ( u Y 



The functions a' v (uS) and a"(u>) are thus identified as the linear suscepti- 
bilities of the change in bond fraction Sv to a small sinusoidal driving force. 
Formally, they have the shape of the susceptibilities of a Kelvin- Voigt material. 
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Note, however, that at the current stage, no mechanical interpretation can be 
made and that the mechanical response of a material vulnerable to internal bond 
breakin g wi ll n ot au tomatically resemble that of a Kelvin- Voigt material (see 
sections 3.1 and 3.2 1. Note further that the frequency scale is set by k- (1 + e u ) 
and not by the equilibrium off rate k alone, as one naively might have expected. 
As an aside, we remark that for our present purpose there is no distinction 
between the unbinding and the partial unfolding^ of the crosslinking protein, 
if the unfolded structure is sufficiently long and floppy, but that there is recent 
evidence from single-molecule experiments 47 that bond-breaking could be the 
more common mechanism for cytoskeletal networks. 



2.2 Power-law fluid 

As a qualitative, analytically tractable model for a transiently crosslinked biopoly- 
mer network, we consider a model power-law fluid that displays weak power-law 
rheologyEl over the frequency range of interest. To enter bond-breaking, we as- 
sume that the linear susceptibility of the power-law fluid is modulated by the 
inverse bond fraction, i.e. 

a' s {u) = G -VH-^- 5 (8) 

and 

atj,(w) = tan (*|) Gq 1 i/(w)- 1 W -', (9) 
with an unspecified prefactor Gq 1 . 




Figure 1: Sketch of a transiently crosslinked polymer network. In thermal equi- 
librium, some polymer-polymer contacts are crosslinked (bright circle segments), 
while other crosslinkers are unbound (dark circle segments). For simplicity, we 
use a reduced description where direct polymer-polymer contacts react by ef- 
fective first-oder kinetics. 



2.3 The inelastic Glassy Wormlike Chain 

The glassy wormlike chain (GWLC) model is a phenomenological extension of 
the wormlike chain (WLC) model, which is the standard model for individual 
scmiflcxible polymers such as DNA or F-actin. For an introduction to the 
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GWLC, see Ref. [THl A detailed description of the inelastic GWLC is given in 
Ref. [35] 

In the weakly bending rod approximation, i.e. if thermal forces do not bend 
the polymer too strongly the WLC can be solved analytically by a mode de- 
composition ansatz for the transverse deflections from the straight ground state. 
The individual modes with index n decay exponentially with a time constant 
r„. The GWLC rests on the phenomenological assumption that the relaxation 
of modes with wavelength A„ larger than the mean bond distance A is slowed 
down, because X n / A— 1 free energy barriers have to be overcome for the mode to 
relax. Formally, this is expressed by a stretching of the mode relaxation times, 

~ _ J T n A„ < A , , 

Tn \_ r n cxp [£(A„/A - 1)] A„ > A " (W) 

In the limit of large free-energy barriers £ , the viscoelastic response of the 
GWLC resembles that of a power-law fluid with exponent^ 3£ _1 . 

In the equilibrium GWLC model, the mean distance A between polymer- 
polymer junctions is taken to be a constant A , unaffected by the equilibrium 
bond fluctuations. The inelastic GWLC keeps track of forced inelastic bond 
breaking in the network by dynamically updating A if the mean fraction v of 
closed bonds in the network changes: A(t) = Ao/i>(t). In the simplest version 
of the inelastic GWLC, the dynamics of the bond network is approximated by 
an effective first-order kinetics and given by the model presented in section |2~Tj 
See figure [T] for an artistic sketch of the physical picture. 

3 Results 

3.1 Apparent linear response of a transiently crosslinked 
power-law fluid 

To gain some analytical insights, we examine the influence of bond breaking on 
the linear response of a schematic power-law fluid, as outlined in section |2.2| 
equations (|8|-(|9j. To this end, we first assume that both the bonds as well 
the material in the absence of bond breaking exhibit a linear response. For 
the sinusoidal driving force in equation Q, the displacement response can be 
written as 

X (t) = G ° 1/V ' 



vq - f ■ sin(w i) - «"( w o) cos(w i)) 

x (sin(wot) — tan(<57r/2) cos(wo^)) , (11) 

where v Q is the steady-state bond fraction in the absence of driving. If we 
expand the right-hand side to first oder in /, we obtain 

x(t) « — (sinfwo^) — tan((57r/2) cos(wot)) , 

i.e. the bond breaking does not contribute to the linear response. 

This suggests that the interpretation for the peaked dissipative patterns ob- 
served in Refs. I31H33I in terms of bond-breaking is not compatible with linear 
response. To see how nonlinear effects can enter the apparent linear response, 
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we note that it is common practice to extract rheological spectra from a weakly 
nonlinear response by Fourier-transforming the response and extracting the con- 
tribution in resonance with the driving to reconstruct a fictitious sinusoidal re- 
sponse with well-defined amplitude and phase shift (see figure [2] for a somewhat 
exaggerated example). As nonlinear processes generally excite higher harmon- 
ics but also change the resonance response, this practice is bound to include a 
nonlinear bycatch. 



a b 




Figure 2: Weakly nonlinear response of the inelastic power-law fluid (red solid 
line) and the corresponding apparent linear response (blue dashed line), recon- 
structed as described in the main text, displayed over driving force (panel a) 
and time (panel b). Parameters are fAxi, — 0.35, 6 — 0.05, lu ~ k. Nonlinear 
effects in the rheology can already be seen for much smaller amplitudes. 

In the following, we determine the nonlinear contributions to the apparent 
linear susceptibilities by calculating the nonlinear response and isolating the 
contributions in resonance with the driving frequency. In order to do this in 
a formally fully consistent way, we would need to include higher-order terms 
for the bond response and for the material response in the absence of bond 
breaking. However, for our present purpose, which is to capture the essential 
physics of the model rather than to deal with some circumstantial details, this 
is in fact not necessary. We therefore assume favorable conditions, namely that 
both the bonds and the material in the absence of bond breaking react linearly, 
and that the only relevant source of nonlinearity lies in their mutual interaction. 
Deviations from our results due to an actually nonlinear bond response and due 
to potential mechanical nonlinearities in the employed viscoleastic model (here 
the power-law fluid) can be evaluated numerically if need arises. Their main 
effect will typically be a mere parameter renormalization. 



With these remarks in mind, we can use equation (111 and expand it to 



higher order in /. Up to third order and using the abbreviation d — tan(<$7r/2), 



6 



the expansion reads 



t 

x(t) « oJq — (sin(cjoi) — dcos(uiot)) 

GqVo 

sffY 

+ G 1 uj a — {sin(wot) — dcos(wot)} {a' v sin(wot) — a" cos(woi)} 

V°J 

- ( f\ 3 

+ G Q 1 uj q d [ — {sm(uj t) — dcos(wo*)} {a' v sin(woi) — a" cos(woi)} , 

(12) 

with the equilibrium bond fraction v§ = (1 + e^ 11 )^ 1 . By rewriting the trigono- 
metric functions one can show that the terms of order / 2 only induce a response 
at the first harmonic at 2wo and a constant offset. Therefore, the first con- 
tribution to the apparent linear response comes from the third-order terms. 
Dropping all terms that do not contribute to the response at coq, the apparent 
linear response becomes 

$ 

ps ojQ i '— (sin(w i) — d cos(woi)) 

Gqv 



'° ~° 4 X 

{ [3a'„(uj ) 2 - a"(w ) 2 - 2da^(w )a"(w )] sin(uj t) 

- [<KM 2 + 3^'(w ) 2 + 2a' u (uj )a'l(ui )] cos(w t)}. (13) 

The apparent linear susceptibility belonging to this response has the real and 
imaginary parts 

a^' M = gL-[l-UL\ [«> ) 2 + 2da> )a>o) - 3a> ) 2 



(14) 



and 



a (3) "M = \d+] (£) [3d<M 2 + 2a> KM +da>o)' 

(15) 

respectively (figure [3]) . Obviously, we can isolate the nonlinear bond-breaking 
contribution by 

, , °' 3) M-° (0) W nfi , 

aW (w ) 

where a is either a' or a" and a^°^(o;o) stands for the linear susceptibilities of 
the material without bond breaking. 



Equation ( 16 ) can be discussed analytically. By inspection of equations ( 14 
16| we can see that exhibits the same scaling behavior as equations (16 
7jT In particular, the frequency scale is also set by A:(l + e u ). However, due 
to the nonlinear combination of a v with as in equation (16 1, the peak of the 
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Figure 3: Nonlinear (symbols) and linearized (lines) susceptibilities (a) and 



moduli (b) for the inelastic a power-law fluid defined in section 2.2 Due to the 
bond breaking, a shoulder develops in the in-phase response (blue), and a peak 
is observed in the out-of-phase response (purple). The parameters are S = 0.01, 
k = 1.4 x 10~ 2 , U = 1, and fAx b = 0.35. 

dissipative nonlinear combined response, (u>), will in general not be located 
at to = k(l + e u ). A simple calculation (see appendix) reveals that the position 
of the peak is at 

u**±=k-{l + e u ). (17) 

This approximate nonlinear formula gives us an analytical prediction (with- 
out any free parameter) of the position of the peak in the spectrum in depen- 
dence of the bond parameters. To assess the quality of the analytical prediction 
for the peak position, we numerically evaluate equations (JlJ-Q (using an Eulcr 
scheme) together with equation ([9| (an example is provided in figure [3| and 
extract the peak frequency for various values of either k or U, keeping all other 
parameters constant. The results are presented in figure |4j We find that the 
functional dependence is predicted correctly, but that the theoretical predic- 
tions are off by a constant factor of about 0.6 (figure [4J3) . This is so because 
we assumed the bonds to always respond linearly and neglected the second- 
order contribution to the bond response, which actually also contributes to the 
third order of the displacement response. If one is interested in quantitative pa- 
rameter estimations, it is straightforward to take into account the second-order 
contribution. However, the more extended formula is rather cumbersome. In 
practical applications, additional nonlinear corrections are more conveniently 
taken into account by a phenomenological correcting factor. Multiplying e u on 



the right-hand side of equation (17) by 0.64 



uj* w -=k ■ (1 + 0.646^), (18) 
v3 

leads to a good quantitative agreement with the numerical results for a broad 
range of k and e u , (figure [Zj dashed lines), independently of the driving force 
amplitude / for fAxb < 0.18 (see supplementary figure 1). 
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Figure 4: Position of the peak in the imaginary part of the nonlinear sus- 
ceptibility. Numerical evaluations for the power-law fluid (blue squares) com- 
pared to numerical evaluations for the inelastic GWLC (red circles, see sec- 
tion 



3.2) and the analytical prediction, equation (TT7J) (solid lines) and the semi 



phenomenological correction, equation ( pL8[ ) (dashed lines, see main text), for 
various values of the equilibrium off rate k (panel a, 10~ 8 < k < 0.1) and the 
equilibrium relative binding affinity e u (panel b, < e u < 64). The driving 
amplitude was fAxb = 0.12, other parameters as in figure [3] For the numerical 
Euler scheme, the time step was chosen as (27r) _1 a;oA< = 5 x 10 -5 . 



The amplitude of the absorption peak turns out to be much more sensitive to 
the intrinsic nonlinearity of the bond dynamics. Deriving a simple approximate 
formula as for the peak position is therefore much less useful. Nevertheless, 
we can predict the approximate dependence on some key parameters from the 
more complicated, fully consistent expressions. Here we only mention that the 
amplitude is independent of the equilibrium off rate k, depends on the relative 
affinity as (1 + e 17 ) -1 and depends quadratically on the driving force amplitude 
/ (data not shown). We further note that it is easy to show from the analytic 
equations that asymptotically, for low frequencies, both the in-phase and the 
out-of-phase part of the susceptibility are merely shifted from the linear-response 
to higher values, where the shift is stronger for the in-phase susceptibility (see 
figure [3J. The interpretation is that the network softens, while the dissipation 
increases due to forced bond-breaking. Note that the high-frequency response is 
not affected by bond breaking, while the low-frequency power-law asymptotics 
violates the Kramers-Kronig relations, testifying the nonlinear character of the 
modulus. 

3.2 Apparent linear response in the inelastic GWLC model 

The inelastic power-law fluid analyzed in the preceding section is to be un- 
derstood as a schematic model for weakly crosslinked biopolymcr networks. A 
more realistic but still minimalistic description is given by the recently proposed 
model of an inelastic glassy wormlike chain^ that builds on the linear modu- 
lus of the phenomcnologically highly successful glassy wormlike chain (GWLC) 



9 



mo( j e llU!6i48| an( j cojnbineg it w jth a bond kinetics scheme as outlined in sec- 
tion |2.3| Due to the slightly more complex structure of the GWLC modulus 
compared to the power-law fluid, we resort to numerical evaluations to explic- 
itly work out the response. We can again compare these numerical results to 
the analytical predictions obtained above to better assess the plausible range of 
validity of our analytical expressions in applications to biopolymer networks. 

As established above for the inelastic power-law fluid, bond breaking matters 
only for the nonlinear response. We therefore again operationally define a set 
of apparent linear moduli by isolating the resonance with the driving frequency 
in the Fourier spectrum of the response. Evaluating the nonlinear response 
of the inelastic GWLC as described in Ref. [35] and calculating the apparent 
linear moduli for oscillatory stress control, we find a qualitative agreement with 
measurements on crosslinked actin networks^ * 33 * 49 !. In particular, we find again 
a maximum in the loss modulus and a shoulder in the storage modulus as 
illustrated in figure [5j 




0.1 1 10 100 1000 



aj/k(l+e u ) 



Figure 5: Apparent linear storage (blue) and loss (purple) modulus versus fre- 
quency (arbitrary units). The lines represent the strictly linear viscoelastic 
storage and loss modulus in the absence of bond breaking. The parameters are 
k = 1.4 x 10~ 2 , U = -1, and fAx b = 0.15. 

It is apparent from figure [5j that the frequency dependence of the linear 
moduli superimposes with the bond effects in a first approximation. There- 



fore, we can easily isolate the nonlinear part of the response (see section 3.1 1. 
To numerically extract the parameter dependence, we evaluate the nonlinear 
frequency-dependent response of the inelastic GWLC at various values of ei- 
ther k or U, keeping all other parameters constant (figures Hand [6}. There 



are essentially two routes to change k (see sections 2.1 and 2.3 1, one assuming 
that changing k also affects the linear modulus of the GWLC (by changing the 
barrier height £) and one that leaves the linear modulus unaffected (assuming 
that only the microscopic attempt rate t^ 1 or the bond stiffness changes). Nu- 
merical evaluations show that both routes lead to the same results for the peak 
positions within numerical tolerance (see supplementary figure 2). In contrast, 
changing U always influences the linear modulus, because U affects equilibrium 
properties such as the steady-state bond fraction. 

The numerical evaluation of the inelastic GWLC in figure [6] clearly reveals 
the dependence of the peak position on the parameters. In the case of a strongly 



frequency-dependent or very large modulus (figkjp), the dissipation due to bond- 
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Figure 6: Apparent linear rheology of an inelastic GWLC. Frequency-dependent 
loss modulus for various values of k (panel a, 10~ 8 < k < 0.1, increasing from 
left to right) and for various values of the relative binding affinity e u (panel b, 
< e u < 10, increasing from bottom to top), all other parameters as in figure [5] 



breaking may be insufficient to become clearly discernible in the rhcological 
spectra, but it can always be uncovered by the subtraction scheme given in 



equation ( 16 ). 

The actual position of the peak in dependence on k (figure [4] a) and e u 
(figure |4]b) obeys the prediction, equation ( 18 1, from the analytical calculations 
in section 3.1 very well (note that we do not have any free parameters to adjust). 
Most strikingly, the numerical values for uj* for the inelastic GWLC agree within 
the numerical error with the numerical values for the inelastic power-law fluid, 
indicating that the dissipative pattern is reasonably robust. Its characteristic 
properties are dominated by the generic bond breaking mechanism and relatively 
insensitive to details of the rheology of the uncrosslinked polymer solution. 

3.3 Exemplary application to actin/ct-actinin networks 



Our analytical prediction for the position of the peak in equation ( 18 1 is sufficient 
to extract valuable kinetic information from bulk rheological measurements. 
While, in principle, one could also exploit the information encoded in the peak 
amplitude, this would require a more sophisticated fitting procedure taking into 
account detailed experimental information on the nonlinearity of the response. 
As this information in not available for pertinent published data, we restrict the 
following analysis to an interpretation of the peak position, which depends on 



the equilibrium off rate and the relative affinity, with equation (18 1. To exem- 
plify the procedure, we use data for an in vitro reconstituted actin/a-actinin 
network^, where a peak in the apparent linear loss modulus was observed. 

Assuming power-law rheology to be a valid model for the (strictly) linear 
response, we first use the high-frequency data for g' to estimate GqVq and 5 to 
4.7 and 0.11, respectively. Next, we use equation (|9| to calculate the linear loss 



modulus . Then, we use equation ( 16 1 to isolate the peak, which turns out 



to be located at w* w 2s 1 . Using an affinity value from the biochemical litera- 
ture^, K — 2.5 x 10 6 M _1 , together with actin and a-actinin concentrations of 
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ca = 23.8/xM and c a = 0.238, respectively^, we obtain e u — K(ca +c a ) ps 60. 

From equation (18l, we find an equilibrium off rate of k ~ 0.09s _1 , which 
agrees with the value obtained for a-actinin/actin in mechanical single- molecule 
experiments^ within the experimental uncertainty. Since our analysis is based 
on a simple (Bell-type^) bond-breaking model the numerical agreement should 
not be over-interpreted, but is certainly reassuring. Note that the rate obtained 
via either route is a factor of five smaller than values obtained biochemically^. 
It has been argued that this discrepancy could be due to the restriction of the 
reaction coordinate by directional loads in single-molecule experiments^. How- 
ever, we noted an interesting fact that may or may not be a systematic feature. 
If we follow experimental hints from actin/HMM solutions that oj* is approxi- 
mately independent of actin- and crosslinker concentration^ and consider the 
concentration values of ca = 2/iM and c a = 3fiM used in the biochemical bulk 
measurement^, we obtain k ps 0.4s -1 , which is exactly the value reported. This 
would suggest that the equilibrium off rate depends on either the actin- or the 
crosslink concentration. Another possible explanation could be that different 
types of a-actinin have been used or that biochemical and biomechanical rate 
measurements systematically differ by a factor, which cancels out when taking 
the ratio of the rates to calculate the affinity. This problem certainly deserves 
a dedicated experimental investigation testing the various hypotheses and pos- 
sibly involving more elaborate models of bond breaking! 43 * 44 !, which is beyond 
the scope of the present contribution. 



4 Conclusions 

We have put forward a simple mechanistic model for the weakly nonlinear rhe- 
ology of reversibly crosslinked biopolymer networks. We verified our analytical 
predictions against numerical results from two generic models for the viscoelas- 
tic response of biopolymer solutions and networks and found good agreement. 
Thereby, we established the characteristic dissipative pattern due to bond break- 
ing as a robust and reproducible rheological signature of transiently crosslinked 
biopolymer solutions with well-defined crossliker kinetics. The model seems 
suitable to explain some puzzling features in the frequency-dependent rheology 
of crosslinked actin networks in terms of the underlying physics and demonstrate 
the non-Maxwellian nature of what could easily be mistaken for a classical ter- 
minal relaxation pattern in a limited frequency interval. 

There are two main conclusions to be drawn from the results presented above. 
First, the non-trivial rheology of transiently crosslinked biopolymer networks as 
observed in Refs. 1321331491 can also be found in a conceptually very simple 
model of an inelastic biopolymer model, the inelastic GWLC. This supports the 
expectation that the inelastic GWLC is indeed a suitable minimal model for 
weakly crosslinked biopolymer networks. Second, as the frequency-dependence 
of the maximum in the loss modulus is very well predicted by our perturbation 
theory, it seems plausible that the mechanism put forward in section |3.1| is 
suitable for deducing information on the crosslinker kinetics from the rheological 
spectra. In doing so, one implicitly acknowledges the intrinsically nonlinear 
character of the absorption peaks revealed by our analysis. 

Compared to single-molecule assays, which provide an alternative method to 
probe crosslinker properties^, rheology has the advantage of being an ensemble 
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method. A potential caveat of the outlined method is that the procedure to 
extract the kinetic information from the modulus relies on the knowledge of 



the linear modulus in the absence of bond breaking (see section 3.1 1. This 
information may sometimes be experimentally inaccessible, e.g. if true linear 
response conditions are technically difficult to realize. However, we discussed 
the nonlinear rheology for two generic models for the linear viscoelastic response 
of biopolymer solutions and networks, and we anticipate that for most practical 
applications at least one of them will provide a reasonable approximation. 

The results presented in this paper thus help to turn a widespread specula- 
tion about the physical mechanism providing biopolymer networks (and possibly 
also living cells) with their extraordinary mechanical properties into a concise 
quantitative prediction amenable to detailed experimental investigation. 
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A Peak position in the nonlinear response 



b = {{Ax b 



Using the abbreviations x = + e u )k, h(x) = c^j b d , 

Ax u )/(4cosh 2 ([//2)) • f /vq) 2 , and d = 7r/2tan(<5), equation (16) becomes 



h{x) 



3bd 



bd 



2(l + x 2 ) 2 4 (l + x 2 ) 



1 + x 2 

The peak position x* is determined by the condition 

d 



2\2 ' 



dx 



h(x) 



= 0. 



(19) 



(20) 



After some simple algebra, equations (19 1 and (20 1 yield the condition 

™*2 



„*3 



_ I * 

d 3 X 



1 

M 



= 0. 



(21) 



For d > 0, this equation has only one positive real solution, x* — l/v3- 
Noting that x* = uj*/(1 + e u )k, we obtain equation ( IT I for the frequency of 
the peak. 
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